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ABSTRACT 



We study the contribution of the kinematic Sunyaev-Zel'dovich (kSZ) effect, gen- 
erated by the warm- hot intergalactic medium (WHIM), to the cosmic microwave 
background (CMB) temperature anisotropies in the Five- Year Wilkinson Microwave 
Anisotropy Probe (WMAP) data. We explore the concordance ACDM cosmological 
model, with and without this kSZ contribution, using a Markov chain Monte Carlo 
algorithm. Our model requires a single extra parameter to describe this new compo- 
nent. Our results show that the inclusion of the kSZ signal improves the fit to the data 
without significantly altering the best-fit cosmological parameters except Q^h 2 . The 
improvement is localized at the I > 500 multipoles. For the best-fit model, this extra 
component peaks at I ~ 450 with an amplitude of 129 uK 2 , and represents 3.1% of 
the total power measured by the Wilkinson Microwave Anisotropy Probe. Neverthe- 
less, at the 2a level a null kSZ contribution is still compatible with the data. Part 
of the detected signal could arise from unmasked point sources and/or Poissonianly 
distributed foreground residuals. A statistically more significant detection requires the 
wider frequency coverage and angular resolution of the forthcoming Planck mission. 

Subject headings: cosmic microwave background. Cosmology: theory. Cosmology: 
observations 



1. Introduction 



Baryons represent a small fraction of the total mass-energy budget of the Universe and do 
not play a predominant role in its evolution. They are the only matter component that has been 
identified directly. The baryon fraction of the Universe has been determined at differen t redshifts 
through a variety of methods: Q h h 2 = 0.020 ± 0.002 (jBurles. Nollett. & Turner] boOll ) from Big 



Institute) de Astrofi'sica de Canarias, Via Lactea, s/n. 38200 La Laguna, Tenerife, Spain; email:rgs@ll. iac.es 
2 Astrophysics Group, Cavendish Laboratory, University of Cambridge CB3 OHE, UK 
3 Ffsica Teorica, Universidad de Salamanca, 37008 Salamanca, Spain; email: atrio@usal.es 

4 Astrophysikalisches Institut Potsdam. D-14482 Potsdam, Germany; email: jpmuecket@aip.de; email: jk- 
lar@aip.de 



- 2 - 



Bang nucleosynthe sis (BBN), fl^h 2 > .021 (|Rauch et al.lll997l ) from the Lya forest and Q^h 2 = 
0.02273 ± 0.00062 (jDunklev et al.ll2009l ) from the cosmic microwave background (CMB) primary 
anisotropies. The baryon fraction m easured from the well-observed components at z = is £l\>h 2 = 
0.010 ± 0.003 (jFukugita et al.lll998l ). indicating that half of the baryons in the local Universe are 
still undetected. 



19951 : iHernauist et al.lll996l : Katz et al.lll996l : iTheuns et al.lll998l : 



Petitiean et al. 


1995: 


Zhang et al. 


Dave et al. 


1999) have shown that 



the intergalactic gas has evolved from the initial density perturbations into a complex network of 
mildly nonlinear filaments in the redshift interval < z < 6. With cosmic evolution, a significant 
fraction of the gas collapses into bound objects; baryons in the intergalactic medium (IGM) are 
in filaments containing Lya systems with low HI column densities and, at low redshifts, shock- 



confined gas with tem peratures T e ~ 0.01-1 keV and overdensities 5b ~ 10 — 50 (jDave et al 



2001 



Cen fc Ostrikerl 12006). An important fraction of the missing baryons could be located in this web 
of shock-heated filaments, called "warm/hot intergalactic medium" (WHIM). Observational efforts 
to dete ct this "missing baryon " component range from looking at its emission in the soft X-ray 
bands (IZappacosta et al.ll2005l ). ultraviolet absorption lines in the spectra of more distant sources 
(jNicastro et al.l 120051 ) or Sunyaev-Zel'dovich (SZ) contributions in the direction of superclusters 
of galaxies (Genova-Santos et al. 2008; see Prochaska & Tumlinson 2008 for a review). Indirect 
searches of the WHIM using the SZ effec t have been inconclusive. Th e SZ imprint due to galaxy 
clusters in WMAP data is well measured (jAtrio-Barandela et al.ll2008br). but when this com ponent 
is removed no signal associated to the WHIM remains (jHernandez-Monteagudo et al.1 12004 ) . 



Atrio-Barandela &: Miicketl (|2006l ) developed a formalism to account for the contribution of the 
IGM/WHIM to the CMB anisotropies via the thermal SZ (tSZ, Sunyaev & Zeldovich 1972) effect. 
The main assumption was that missing baryons were distributed as a diffuse gas phase outside 
bound objects. Its filamentary structure was assumed to be described by a log-normal distribution 
function, which accu rately models a mild ly non-linear density field when the velocity field remains 
in the linear regime (jColes &; Joneslll99ll ). The predicted power spectrum of the tSZ effect peaks 
at I ~ 2000-4000, with an am plitude similar to the tSZ from galaxy clusters. More recently, 



m 



Atrio-Barandela et al.l (|2008al ) we studied the contribution of the kinematic SZ (kSZ, Sunyaev 
& Zeldovich 1980) effect. We found that the kSZ power spectrum had a maximum at £ ~ 400- 
600. In this article we search for a possible kSZ contribution in the 5 year Wilkinson Microwave 
Anisotropy Probe ( WMAP) data. We use a Markov chain Monte Carlo (MCMC) method to sample 
the parameter space of the concordance ACDM model with and without a kSZ component to 
determine if there is a statistically significant contribution. Briefly, in Sections 2 and 3 we describe 
the model and the numerical implementation of the MCMC, and in Section 4 we present our results 
and summarize our main conclusions. 
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2. The thermal and kinematic SZ effect from the IGM/WHIM. 



The tSZ effect is the weighted average of the electron pressure along the line of sight; the kSZ 
is proportional to the column density of the free electrons along the line of sight h, weighted by the 
radial component of their peculiar velocities: 
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(1) 



In these expressions, T e , n e , v e are the electron temperature, density and peculiar velocity, respec- 
tively, k-Q is Boltzmann constant, <tt Thompson cross section, m e c 2 the electron annihilation energy, 
c is the speed of light and G{v) is the frequency dependence of the tSZ effect. In the Rayleigh- Jeans 
regime, G{y) ~ —2 with less than 20% variation at WMAP frequencies. For an isothermal cluster 
the tSZ to kSZ ratio is: (AT t sz/AT kSZ ) ~ 20G(i/)(T e /10keV)(300kms- 1 / u e)- The temperature 
of the IGM is much lower than in galaxy clusters even in the shock-heated WHIM and the kSZ 
contribution could become comparable to that of the tSZ effect. 



In lAtrio-Barandela Sz Miicketl (j2006l ) and lAtrio-Barandela et al.1 (|2008al ) we computed the tSZ 
and kSZ temperature anisotropies generated by the IGM assuming that baryons are distributed 
as in a log-normal random field. The log-normal distribution was introduced by Coles & Jones 
(1991) as a model for the non-linear distribution of matter in the Universe. The number density 
of electrons n e can be obtained from the baryon distribution assuming ionization equilibrium be- 
tween recombination and photoionization and collisional ionization. In the conditions valid for the 
photoionized IGM the gas is almost completely ionized. The correlation function of the tSZ tem- 
perature anisotropy generated by two filaments located at two different redshifts along two lines of 
sight with an angul ar separation a is dominated by t he spatial variations of the electron pressure 
at nearby locations (lAtrio-Barandela &; Miicketl 12000 ) . Within the small angle approximation: 
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In this expression, A 2 (z) is the variance of the baryon density field, which is related to the dark 
matter power spectrum P(k) by: 



A\z)=D 2 (z) J 



d?k P(k) 
(2^) 3 [1 + x 2 b (z)k 2 ] 2 



(3) 



where D(z) is the growth factor of matter density perturbations and Xf, is the comoving Jeans 
length. Also, 
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2vr 2 J [l+x 2 {z)k 2 ] 2 

where jo is the zeroth order spherical Bessel function. The integration in eq. ([2]) extends to the 
highest redshift zj. 
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Similar arguments can be used to compute the correlation function of the kSZ effect due to 
two filaments located at dista nces I and /' along two lines of sight separated by an angle a (see 



Atrio-Barandela et al.l (|2008al ) for further details): 



2 r i(z f ) r V(z f ) 

C 2 



C(a) = fZ^r / dldl'VB(R(z))D v (z)VB(R(z'))D v (z')(n e (lh,z)n e (l'n',z')) , (5) 



where Vb(R(z)) denotes the mean bulk velocity of a sphere with electron density n e and radius 
R(z), the comoving distance from the observer to a filament at redshift z, D v (z) is the velocity 
linear growth factor and /b is the fraction of the baryons in WHIM filaments. 

The correlation functions given in eqs ([2]) and ([5]) differ in two significant aspects: (a) While 
the tSZ effect is cumulative and several filaments along the line of sight add linearly to the total 
effect, eq. ([5]) requires all filaments to be moving with the same velocity. This is the physical reason 
why only the bulk flow velocity contributes to the effect. If there are several filaments, in a given 
direction, with different velocities the net effect will not be cx (n e T^) but oc (n e V B ) 1/2 . (b) TheIGM 
is usually described by a polytropic equation of state T oc n 7_1 , i.e., high dense regions contribute 
more to the tSZ effect since their temperature is higher. In eq. ([2]) we then integrate all scales out 
to the largest overdensity (5 max that is well described by the log-normal model. Otherwise we would 
be including small contributions from gas in bound objects and not described by the log-normal 
model. This restriction is not necessary for the kSZ effect. In the log-normal approximation, 
velocities are in the linear regime and they are not correlated with matter overdensities. Then, 
instead of introducing an arbitrary cut-off at overdensity <5 max , we extend our integration to all 
overdensities, i.e., to all baryons, and introduce a factor /b to take into account the actual fraction 
of baryons in the WHIM; we could then, in principle, constrain this fraction directly from the data. 

Eqs ([2]) and ([5]) can be inverted to give respectively the radiation power spectra C^ sz and 
Cf sz at each multipole I. These calculations are computationally expensive since, to get accurate 
results, we use line-of-sight separations of 30" and redshift intervals of Az = 0.005 up to z = 0.1 
and Az = 0.1 up to Zf = 1, where the contribution becomes negligible. The final tSZ power 
spectrum depends on cosmological and physical parameters: the cut-off scale 5 max , the amplitude 
of the matter density fluctuations on spheres of 8/j" 1 Mpc, ag, the mean gas temperature, To, the 
gas polytropic index, 7, and the mean gas temperature at reionization, T m . The kSZ contribution 
depends on the Jeans length x&, <7s and /b- Since the product jT m fixes xj, then the kSZ contribution 
requires less parameters than the tSZ. 

In Figure QJ, we compare CMB power spectrum of the concordance ACDM model with the 
kSZ and tSZ WHIM contributions. The parameters of the concordance model are those of the 
best-fit 5 year WMAP data. The SZ contributions are calculated using: og = 0.77, 7 = 1.3, 
7T m ~ 1.3 x 10 4 K and /b = 0.5. First, the relative amplitude of the thermal and kinematic 
contributions depends on model parameters; the kSZ effect could be larger than the tSZ effect 
since the average temperature of the IGM is rather low. Second, since the IGM is not isothermal, 
the tSZ contribution is dominated by the mildly overdense regions, which subtend smaller angles 
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than the filaments themselves, so its power is shifted to higher I. In Figures [Th>,c we show the 
contribution of different redshift intervals to the total power for each of the WHIM contributions. 
In both cases, the main contribution comes from z = to 0.4. Contributions from higher redshifts 
decrease rapidly. 

The purpose of this article is to search for a IGM/WHIM contribution in the 5 year WMAP 
data, which are sensitive to I < 1000 multipole s. On these angular sc ales, the kSZ contribution 



is largest while the tSZ is larger at £ ~ 3000. Friedm an et al.l (|2009l ) see no evidence of a tSZ 
component at 2000 < I < 3000 and I Sharp et al.l (|2009r ) constrain the tSZ power spectrum at 
I ~ 4000 to be < 149 /iK 2 at 95% confidence level. The expected tSZ contribution at WMAP 
scales is then negligible. For this reason, we shall consider only a kSZ component. This restriction 
simplifies our study since the kSZ power spectrum depends only on three parameters: as, jT m and 
/t>. In Figure [2] we plot kSZ power spectra for different values of as and ^T m with /b = 0.5. A 
regression fit of the variation of the maximum amplitude of the kSZ power spectrum with each of 
these three parameters permits us to write the following scaling relation: 
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While the amplitude depends on the three parameters, the location of the maximum is only 
weakly dependent on jT m . Since /b is a multiplicative factor that affects only the amplitude but 
not the shape or position of the maximum, variations on /b and as cannot be distinguished. The 
kSZ effect is then effectively described by two parameters: one cosmological (as) and one physical 
(jT m ) determining the Jeans length. 

Figure [Tfc also shows that most of the contribution comes from very low redshifts. Similar 
scaling relations hold for the tSZ contribution of clusters of galaxies, indicating that we could 
generate arbitrary large contributions. I n semi-analytical estimates, the abundance of clusters is 



1999 



Molnar &: Birkinshaw 



given by the Press-Schechter formalism (jAtrio-Barandela Miicket 
2000) and a large tSZ effect is obtained by arbitrarily increasing the number of clusters. The 
scaling relation given in eq. ([6]) is limited by the validity of eq. ([5]); if filaments overlap along the 
line of sight, it overpredicts the signal. Since the kSZ effect is just an integral along the line of 
sig ht of all elect rons with a coherent peculiar motion, an order of magnitude estimate of the effect 
is (|Hoganlll992h : 

$b \ ( fb \ ( L \ ( Vb 



(AT) 



kSZ 



' /,K ' 50 J 1 0.5 



(7) 



^30 Mpc J V600 km/s, 

where <5b is the overdensity o f a typical filament and L is the coherence scale of a motion with 
amplitude Vb- As explained in lAtrio-Barandela et al.l (|2008al ). the distribution of the effect is rather 
skewed and 9% of all lines of sight will produce an effect 1.8-6 times larger than the estimate from 
eq. ([7]). In the concordance model, bulk flow velocities of ~ 200 km/s are typical for volumes 
of R ~ 100 Mpc/h radius. Kashlinsky et al. (2008, 2009) reported a bulk flow of amplitude 
600-1000 km/s on a scale of 300 Mpc/h that could give a much larger contribution. A significant 
fraction of the temperature decrement of —230 /uK detected in the intercluster medium of the Corona 
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Borealis Supercluster by iGenova-Santos et al.l (|2005l ) could be due to thermal and kinematic SZ 
contributions. However, that decrement is in a direction almost perpendicular to the bulk flow 
cited above and, because of its orientation with respect to the line of sight, this flow would not 
induce a significant kSZ contribution in Corona Borealis. Only a smaller kSZ contribution could 
exist due to peculiar motions on the scale of the supercluster itself. 



3. Markov chain Monte Carlo parameter estimation 



To explore the parameter space of ACDM mo dels with and without a kSZ contribution, we used 
the April 2008 version of the COSMOMC package (jLewis fc Bridle! |2002| ). This software implements 
an MCMC method that performs parameter estimation using a Bayesian approach. When the kSZ 
contribution is included, precomputed Cf sz are added at ea ch step of the chain to the theoretical 
power spectrum, computed with CAMB (jLewis et al.l [2000j) for each set of ACDM cosmological 
parameters. The model is compa red with the data using the likelihood code supplied by the 
WMAP team (jDunklev et al.ll200g| ). 



We considered the concordance ACDM model, defined by a spatially flat Universe with cold 
dark matter (CDM), baryons and a cosmological constant A. The relative contributions of these 
components are given in units of the critical density: £l c dm, an d = 1 — ^m> where fi B = 
^cdm + is the total matter density. The specific parameters used in the analysis were the 
physical densities fl^h 2 and Q C( i m h 2 where h = Hq/100 km s _1 Mpc -1 is the normalized Hubble 
constant. In order to minimize degeneracies, instead of h we used the angular size of the first 
acoustic pe ak 0, i.e. the ratio of the sound horizon to the angular diameter distance to last 



scattering (IKosowsky et al.1 |2002j). We considered adiabatic initial conditions and assumed an 
instantaneous reionization parameterized by its optical depth to Thomson scattering up to the 
moment of decoupling r. The initial fluctuation spectrum was parameterized as a power law, 

P(k) = A S \JL) , (8) 

where A$ is the amplitude at k c = 0.05 Mpc -1 and n s the spectral index. Since the effect of /b on 
the kSZ power spectrum is indistinguishable from that of ag we fix that parameter to a given value. 
Thus, the kSZ signal is defined by as, which is derived from the previous cosmological parameters, 
and "fT m . In summary, our cosmological model with the kSZ component has seven degrees of free- 
dom and is described by the parameterization 6 = {Qbh 2 , ^cdm^ 2 > 1000, ln(10 10 As),n s , r, log(7T m )]. 

In Table Q] we give the flat priors imposed on each parameter, the initial values and distribution 
widths. In addition, we used a top-hat prior 10 < to < 20 Gy for the age of the U niverse. When 
running the chains we included only the 5 year WMAP data (jHinshaw et al.ll2009l ) and no other 
CMB or cosmological datasets. 

At each step in the chain, a set of initial values for the cosmological parameters and for "fT m 
are generated. The ACDM radiation power spectrum and a$ are computed for the basic model. 
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When introducing the kSZ component into the model, we would have to compute the values of Cf sz 
at each step of the chain. However, as the computation is very demanding, we only calculate kSZ 
power spectra on a 2D eg — jT m grid. We built a grid with a bin width Ads = 0.05 in the interval 
a 8 = [0.60,0.90] and A( 7 T m ) = 1 x 10 4 K in the interval jT m = [1,5] x 10 4 K (plus interleaved 
values at 0.8, 0.9 and 1.3xl0 4 K). At some random points in the parameter space we checked that 
the interpolated and exact spectrum did not differ by more than 7% in the interval £ ~ 100-600. 
After a first run, we resampled our grid around the best-fit model. In the intervals as = [0.70, 0.83] 
and jT m = [1.05, 1.95] x 10 4 K, the bin widths were Aa 8 = 0.01 and A( 7 T m ) = 0.05 x 10 4 K, 
respectively. With this resampling, the differences between interpolated and exact power spectra 
were lower than ~ 2% in the same £ range. The kSZ power spectrum is obtained at each step of 
the chain by logarithmic interpolation of the precomputed power spectra in the four closest nodes 
of this grid. 

Initially, between two consecutive steps in the chain, we added to each parameter a random 
increment drawn from a Gaussian distribution with a standard deviation equal to the value listed 
in Table [1] multiplied by 2.4. Since we found that 7 T m was strongly correlated with Q^h 2 and 
ln(10 10 ^4s), we reran the chains using the eigenvalues of the covariance matrix of the parameters 
(computed with the getdist facility which is part of the COSMOMC package) as distribution widths. 
We ran nine independent chains, with a total number of 150 000 independent samples, when no kSZ 
is included. With kSZ we fixed the baryon fraction in WHIM at five different values /t, =0.3, 0.4, 
0.5, 0.6 and 0.7 (we also tried to constrain a model with /b as a free parameter, but the degeneracies 
did not allow us to obtain reliable results). In each case we ran eight inde pendent chains of similar 



sizes with a total number of samples > 225 000. We used the R statistic (jGelman Rubinlll992l ) 
as a convergence criterion. All our parameters have R well below 1.2: R ~ 1.008 for log( 7 T m ) and 
«1.002 for the other parameters. 



4. Results and discussion 



In Figures [3^1, b,c we plot the mean ID likelihoods for 7 T m , as and the more informative rms 
temperature fluctuation introduced by the kSZ signal defined as: 
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The different lines correspond to different baryon fractions. Since the kSZ parameters affect the 
amplitude of the spectrum most significantly (see Figure [2]) then, for each the maximum ID 
likelihood of 7 T m shifts but the kSZ signal remains roughly constant (see Figure [3b). This high 
degeneracy indicates that 5 year WMAP data is insensitive to the fraction of baryons in the WHIM 
and hereafter we will quote results for /b = 0.5. 

In Figures [3ti,e,f we show the mean 2D likelihoods for pairs of parameters. Solid lines represent 
la and 2a contours, whereas the white dots indicate the positions of the maximum likelihood in 
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the full parameter space. Not unexpectedly, their positions are slightly shifted with respect to the 
2D likelihood maxima. In fact, Markov chains estimate confidence intervals more precisely than 
locate the ma ximum of the likelihood, the bias being larger the greater the model dimensionality 



(jLiddld 12004 ) . Figure [3H shows a noticeable degeneracy between as and log(7T m ). Within the 
la confidence region ^4ksz varies between ~ 20 fjK 2 and ~ 300 /xK 2 . Re gions with high values of 
o"8 and low values of log(7T m ) are strongly ruled out since they overpredict the kSZ contribution. 
At the 2a level absence of a kSZ contribution is compatible with the data. This is also seen in 
Figures EU,e. Models with kSZ prefer a lower value of as (see Figure [3b), since adding kSZ requires 
less primordial CMB power. Figure [3b(f) shows that ^}\,h 2 is (inversely) proportional to (^j^) 1 ' 2 
(log[7T m ]). Physically, kSZ adds power mainly at t ~ 400-600 and the amplitude of the second 
acoustic peak gets reduced to fit the data. As a result, the best-fit value for Q^h 2 increases. 

The parameters that maximize the likelihood are given in Table [21 The la confidence intervals 
have been computed from the ID mean likelihood distri butions. The values for the concordance 



model agree well with those given by the WMAP team (jDunkley et al.ll2009l ). indicating that we 
are correctly exploring the parameter space. The model with kSZ has larger error bars, as we 
are fitting the same data with an extra degree of freedom. The best-fit model has a kSZ power 
spectrum with a maximum amplitude of ^4ksz = 12Q~t}^ M-^- 2 centred at I « 438, corresponding to 
jTm = 1.3^Q2 x 10 4 K. The rms temperature fluctuation is (T 2 ^) 1 ^ 2 = 19Zq /iK, very significant 
compared with the contribution from the primordial CMB, (T^g) 1 / 2 = 110 fiK. None of the 
parameters except tl^h 2 differ by more than la from those of the concordance model. The model 
with kSZ gives a higher value for Q.\>h 2 . This reinforces the slight discrepancy between the values 
obtained from CMB temperature anisotropics and from BBN. As indicated above, adding kSZ 
boosts the value of Q\>h 2 since it reduces the height of the second acoustic peak. The difference in 
as is also notable even if the confidence regions overlap at the la level. Our estimate is closer to the 
values obtained from the number density of gala x y clu sters and the optical or X-ray cluster mass 



functions. For example, IVoevodkin &: Vikhlininl (|2004l ) found as = 0.72 ± 0.04 from the baryon 



mass fraction of a sample of 63 X-ray clusters. After compiling cluster determinations since 2001 



Hetterscheidt et al.l (|2007l ) obtained as = 0.728 ± 0.035, whereas from another compilation of the 



weak lens cosmic shear they found as = 0.847 ± 0.029. 

In Figure Hh. we compare the accuracy of the fitting of the two models to the 5 year WMAP 
data. We plot the best-fit power spectra for the models with and without kSZ after subtracting the 
WMAP band powers. The kSZ model achieves a better fit to the experimental data points in the 
multipole range £ ~ 300-700. The same conclusion can be obtained from Figure Hb, where we plot 
the ratio of the binned x 2 P er ^-band for the model without and with kSZ (note that these values 
were calculated, just for illustrative purposes, by applying the x 2 statistics to the binned best-fit 
theoretical power spectra, and not by the WMAP likelihood code). This ratio is overall > 1 in 
the range i ~ 300-700, even though within this interval there are some particular bins where the 
concordance model produces a slightly better fit. 

Introducing the kSZ component, the x 2 reduces by Ax 2 = —3.3, from 2661.05 to 2657.83. 
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Evaluating the statistical significance of this result requires taking into account the number of 
degrees of freedom, given by the difference between the number of independent data points N and 
of model parameters k. The WMAP likelih ood is computed as a sum of different temperature 
and polarization terms (jDunklev et al.l 120 09). In this computation 968 points correspond to the 
TT power spectrum at t = 33-1000 and 427 to the TE cross-correlation at I = 24-450. For low 
(I < 23) multipoles the likelihood associated with the TT, TE, EE and EB correlations is evaluated 
directly from 1170 pixels of the temperature and polarization maps. In total N = 2565. Adding 
the kSZ contribution raises the number of parameters from k = 6 to k = 7. Note that the x 2 P er 
degree of freedom also decreases when the kSZ component is included, from x^ of = 1.0399 to 1.0390. 
Increasing the number of model parameters always improves the fit to a particular dataset but the 
model loses predictive power. In Bayesian statistics, information criteria can be used to decide 
whether the introduc tion of a new parameter is favored by the data. Examples are the Akaike 
information criterion (lAkaike 1974 ) defined as AIC = 2k + x 2 or the more conservative Bayesian 
information criterion ( Schwarzlll978l ) BIC = x 2 + klnN. Since we are adding a single parameter, we 
obtain AAIC = —1.2, which is marginal evidence in favor of introducing this new parameter while 
ABIC = 4.1 is evidence against it, reflecting the fact that no kSZ component is compatible with the 
data at the 2a level. Even if at present model selection doe s not clearly favors the k SZ contribution, 
it is a model well motivated physically and, as remarked by lLinder Miquell (|2008l ). model selection 
needs to involve physical insight. A statistically more significant measurement would require a wider 
frequency coverage and angular resolution, as will be provided by the forthcoming Planck mission. 

The kSZ contribution represents 3.1% of the power measured in 5 year WMAP data. This con- 
tribution is larger than expected; typical filaments would give rise to ATksz ~ 2-5 (eq. [7]). We 
would need a more accurate (numerical) model to establish which would be the physical conditions 
that give rise to the quoted kSZ signal. Numerical simulations over large cosmological volumes are 
computationally very expensive because they require high resolution over most of the volume. In 
fact, spatial resolution is usually an important limitation. Adaptive mesh refinement (AMR) tech- 
niques are efficient at describing the dyn amics of gas in the high de nse regions but their resolution 
falls sharply in less dense environments (jRefregier &; Tevssierl I2002Q . while smooth particle h ydro- 
dynamics (SPH) codes do not yet resolve scales as small as the Jeans length (Bond et al.ll2005r) . that 



effecti vely dominate the WHIM S Z contribution (lAtrio-Barandela &: Mucketll2006l ). lHallman et al 



(j2007j) used the AMR Enzo code (jO'Shea et al.ll2005l ) to simulate a (512 Mpc/h) 3 volume including 
unbound bas, and found that one-third of the SZ flux in a 100 square-degree region comes from 
objects with masses below 5.0 x 10 3 and filamentary structu r es ma d e up of WHIM gas, a 
conclusion similar to that reached by iHernandez-Monteagudo et al.1 (|2006l ). lHallman et al.1 ((2009) 
focused in the low-density WHIM gas by restricting their analysis of the same AMR simulation to 
regions with temperatures in the range 10 5 — 10 7 K and overdensities 5 < 50, even though their 
low mass halos were not yet resolved gravitationally. They computed the radiation power spectrum 
from that simulation and found a similar shape to that of Figure [Vj although with a much smaller 
amplitude. The hig h amplitude we found c ould be an indirect confirmat ion of large scale peculiar 
motions reported in iKashlinskv et al.l (|2008l ) and iKashlinskv et al.l (J2009J) . If bulk flows of this am- 



- 10 - 



plitude are very common, filaments at higher redshift could also contribute significantly, thereby 
increasing the total signal. Also, we must take into account a possible foreground contamination: 
power spectra with a single maximum are also good models for Poissonianly distributed foreground 
residuals or unresolved point sources. Since the WMAP window function exponentially damps 
power at I > 600, the convolution of a Cg = const, spectrum with this window function peaks at 
(. < 600 and has a shape similar to that of the spectrum de scribed above. Fore ground substraction 



can differ in the range 3-5 /iK depending on the method (jGhosh et al.1 12009 ) and some residuals 
could be present on the data at that level. To distinguish the WHIM kSZ and foreground signals 
we will have to include the frequency dependence, amplitude and location of the maximum for each 
component. In the simplest model with kSZ and one foreground, we would need four parameters 
to model both components so an analysis of all components is unfeasible with WMAP 3 frequency 
bands (Q, V and W). 

To conclude, we have explored the parameter space of the concordance model to show that 
the WHIM kSZ contribution could be as high as 3% the total power of 5 year WMAP data. This 
large amplitude is difficult to account for in the concordance ACDM model, where filaments are 
expected to have < 5 /uK contributions. It could be an indication that large scale flows are rather 
common. We cannot rule out that part of this contribution could be due to unmasked point sources 
and/or foreground residuals. The Planck satellite, with its wider frequency coverage, lower noise 
and different scanning strategy, is well suited for detecting the IGM/WHIM thermal and kinematic 
contributions with much higher statistical significance and to distinguishing this signal from other 
foreground contributions. 

We are thankful to Rafael Rebolo for enlightening discussion and suggestions. This work is 
supported by the Ministerio de Educacion y Ciencia and the "Junta de Castilla y Leon" in Spain 
(FIS2006-05319, programa de financiacion de la actividad investigadora del grupo de Excelencia 
GR-234). 
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Fig. 1. — (a) Intrinsic CMB radiation power spectrum (solid line) and the IGM/WHIM tSZ and 
kSZ (dashed lines) contributions, (b) Contributions of redshift intervals to the total tSZ spectrum. 
From left to right, solid lines correspond to contributions from the following redshift intervals: 
[0,0.005], [0.005,0.05], [0.05,0.1], [0.1,0.2], [0.2,0.3] and [0.3,0.4]. (c) Redshift contribution to the 
kSZ effect. Solid lines give the contribution for z = to 0.1 with spacing Az = 0.01; dotted lines 
give the contributions for z = 0.1 to 0.4 with Az = 0.1. Dashed lines show the total power spectra 
from all the redshift intervals. The SZ contributions are calculated in all cases using as = 0.77, 
7 = 1.3, 7 T m = 1.3 x 10 4 K and f h = 0.5. 




Fig. 2. — Kinematic Sunyaev-Zel'dovich power spectra from WHIM for a fixed as = 0.80 and 
different values of ^T m (solid lines), and for a fixed ^T m = 1.3 x 10 4 K and different values of as 
(dotted lines). The fraction of baryons stored in the form of WHIM has been fixed at /b = 0.5 in 
all cases. 
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Fig. 3. — Top: mean ID likelihoods for the Jeans length 7T m , as and the rms temperature 
anisotropy contribution of the kSZ component. Each line correspond to a different fraction of 
baryons in the form of WHIM. In (b) we also plot the likelihood for the concordance (no kSZ) 
model. Bottom: mean 2D likelihoods for three different combinations of parameters. Solid lines 
depict the la and 2a confidence regions, and the thick dots indicate the positions of the likelihood 
maxima found in the full parameter space. Dotted curves in (d) represent levels with the same kSZ 
amplitude. 
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Basic 


Parameter 


Starting 


Distribution 


parameter 


limits 


points 


width 


Q h h 2 


(0.005, 0.1) 


0.0223 


0.001 


O h 2 
J 'cdm' 1 


(0.01, 0.99) 


0.105 


0.01 


1000 


(0.5, 10) 


1.04 


0.002 


ln(10 10 A s ) 


(2.7, 4.0) 


3.0 


0.01 


n s 


(0.5, 1.5) 


0.95 


0.01 


T 


(0.01, 0.8) 


0.09 


0.03 


log{jT m ) 


(3.778, 4.740) 


4.096 


0.006 



Table 1: For each independent parameter, we indicate its range of variation (top-hat prior), the 
initital values and the estimated distribution widths initially used by the chain to vary each pa- 
rameter. 
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Fig. 4. — Top: best-fit power spectra for the models with (solid line) and without (dotted line) 
kSZ, after subtracting the WMAP experimental band powers. The vertical lines represent the la 
error bars of the 5 year WMAP data. Bottom: ratio of the \ 2 m ^-bands of width A£ = 15 for a 
model without and with a kSZ component. 
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Parameter 
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CMB and kSZ 


n h h 2 

O h 2 

1006* 

ln(10 10 ^ s ) 


n n99^+ - 0007 
u.u^zo_ 0006 

n i nsfi+ - 0081 
0.1U8b_ 00060 

1 040+ - 004 

1.U4U_q QQ3 

3.06 ± 0.05 


n234+ - 0020 
in96 +0 - 0068 

u.iuyo_ 0085 
1 n42+ - 005 

1.U4Z_q 004 

o D9+0-05 
°- uz -0.08 


n s 


0.96 ± 0.02 


0.95 ± 0.02 


T 

7 T m (10 4 K) 


o.o88±8:SS 


0.0871K 

1 01 +0.41 

i -° 1 -0.23 


-Zre 

Ho 

(^ k 2 sz) 1/2 M 
^ksz (^K 2 ) 


0.74518;^ 
u.zoo_ 027 

7OQ+0.042 

u./yo_ 037 
71.63±|;f| 


n 7co+0.043 
u. <oz_ Q 034 

0.2481^ 

n 771 +0.039 
u - ' ' i -0.069 

10.141};" 
73.18l|^ 
19.5l|? 
128.7™° 


x 2 


2661.05 


2657.83 



Table 2: Best-fit parameters for the model without and with a kSZ contribution. These results 
correspond to a fraction of baryons in the form of WHIM of /b = 0.5. The upper and lower (below 
the line) sets correspond to independent and derived parameters, respectively. The central values 
have been derived from the sample with the minimum \ 2 i n the chains (also shown in the table), 
whereas the la confidence limits were derived from the 0.159 and 0.841 points of the cumulative 
probability distribution given by the ID mean likelihoods. 



